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Abstract 

In this paper, we study the Wong-Zakai approximation of the solution to the stochastic 
differential equation with reflecting boundary condition on a domain D in a Euclidean space. 
We prove the L p convergence of the approximation in C([0, T] — > D) under some general 
conditions on D. 

(N ! 1 Introduction 

> 

Stochastic differential equations(SDEs) are defined as stochastic integral equations. The def- 
inition of the stochastic integrals is based on martingale theory although there are pathwise 
approaches to this problems via rough path theory recently. A simple relation between SDE and 
usual ordinal differential equation(=ODE) were found by Wong and Zakai [12]. That is, they 
consider Stratonovich SDE and corresponding ODE which is obtained by replacing the Brownian 
motion by the piecewise linear approximation and prove that the solution of the ODE converges 
to the solution of the Stratonovich SDE almost surely in the topology of uniform convergence 
when the approximation becomes finer. More general approximations of paths are found, e.g., in 
[I]. When we consider SDE on a domain D in M. d , we need to consider boundary conditions. In 
this paper, we study Wong-Zakai approximations of solutions to SDE with reflecting boundary 
conditions on D and prove the LP convergence of them to the solution in C([0,T] — > D). This is 
not a first study of Wong-Zakai approximation of reflecting SDE. Doss and Priouret [2] proved 
the uniform convergence of the Wong-Zakai approximations in probability in the case where dD 
is sufficiently smooth. Also Pettersson [6] proved the uniform convergence in the case where D 
is bounded convex and the diffusion coefficient is a constant matrix. More recently, Evans and 
Stroock [3j proved the weak convergence of the law of the Wong-Zakai approximations for more 
general domains which satisfy conditions (A), (B) and (C) and the admissibility condition. We 
explain these conditions in Section 2. Our results improve their weak convergence to LP con- 
vergence in C([0, T] — > D). We note that there are studies of Euler approximations of reflecting 
SDE. We refer them to [EH IS]- 

The paper is organized as follows. In Section 2, we state our main theorem (Theorem 12. 9p . 
First, we recall the basic results on the Skorohod problems and the existence and uniqueness 
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of the strong solutions of reflecting SDE based on Lions and Sznitman [5] and Saisho [5J. In 
particular, we explain the conditions (A), (B) and (C) on domains following [8j. In Section 
3, we prove L p convergence of Euler approximations under conditions (A), (B) and (C) in 
C([0, T] — t- D). In Section 4, we prove our main theorem. 

2 Preliminary and main theorem 

Let D be a connected domain in R . In this paper, we do not assume the boundedness of the 
boundary of D or D itself. We define the set M x of inward normal vectors at the boundary point 
x G dD by 

M x = u r>Q Af Xi r (2.1) 
M x ,r = {n G M d I |n| = 1, B(x - rn, r) n £> = 0} , (2.2) 

where B(z, r) = {y G l d | |y — z| < r}, z G M d , r > 0. 
Let us recall conditions (A), (B), (C) following [8]. 

Definition 2.1. (1) Condition (A) (uniform exterior sphere condition) . There exists a constant 
ro > such that 

= .A4.ro / ^ for an y » G <9D. (2.3) 

(2) Condition (B). There exist constants 5 > and (3 > 1 satisfying: 
for any a; G 3D there exists a unit vector Z x such that 

(4,n) > i for any n G U 2/ejB ( a . j(5)n9D A/' l/ . (2.4) 

(3) Condition (C). There exists a function on R d and a positive constant 7 such that for 
any x G 3D, y G D, n G A/" x it holds that 

(y - x, n) + - ((D/)(x), n) \y - x\ 2 > 0. (2.5) 

7 

In this paper, the function space denotes a set of C k functions such that all their deriva- 
tives and themselves are bounded. Note that if D is a convex domain, the condition (C) holds for 
/ = 0. The admissibility condition on D in [5| is the property that D can be approximated by 
domains with smooth boundary in a certain sense. In this paper, we do not use such a property 
and we refer it to [5]. Here we explain what Skorohod problem is. Let w = w(t) (0 < t < T) 
be a continuous path on M. d with w(0) G D. The pair of paths (£, <p) on M. d is a solution of a 
Skorohod problem associated with w if the following properties hold. 

(i) £ = (0 < t < T) is a continuous path in D with £(0) = w(0). 

(ii) It holds that f (i) = 10 (i) + for all < t < T. 
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(iii) cf) = cj)(t) (0 < t < T) is a continuous bounded variation path on W 1 such that 0(0) = 
and 

0(t)= r n ( S )dii0ii [OiS] (2.6) 







[o,t] = / l0D(£(s))d|MI[o, s ]- (2-7) 

JO 

where n(i) G 7V e(t) if £(t) G 

In the above, 1 1 1| ro,t] stands for the total variation norm of (f>. See (|2,10|) . 

The existence and uniqueness of solutions were proved by Tanaka [11] for the convex domain 
with additional assumptions. Lions and Sznitman proved the existence and uniqueness under 
conditions (A), (B), (C) and admissibility condition. The following result is proved by Saisho[8]. 

Theorem 2.2. Assume conditions (A) and (B). Then there exists a unique solution to the 
Skorohod problem for any continuous path w. Moreover the mapping V : w t— > £ is continuous 
in the uniform convergence topology. 

Doss and Priouret [2] proved the convergence of Wong-Zakai approximation. They used the 
Lipschitz continuity of the Skorohod map r : w >— > £ in the half space case. Under conditions (A) 
and (B), it is proved that T is l/2-H61der continuous map in the uniform convergence topology. 
See [5]. If T is Lipschitz continuous, Doss and Priouret's approach may be applicable. We 
may use L(w) = T(w) — w which corresponds to the local time at the boundary 3D. 

The bounded variation norm of <p can be controlled by the supremum norm of w and the 
modulus of continuity. Such an estimate is proved by Tanaka [IT] in the case of convex domains. 
Similar estimates are obtained by Saisho [8] without assumptions of the convexity. For our 
purpose, we need quantitative version of Saisho's estimate. To this end, we introduce norms of 
continuous paths. Let < 6 < 1 and define the Holder norm of the continuous path w G W by 

\w(v)-w(u)\ 

IHkM,*/2= /up ■ m ■ (2-8) 

s<u<v<t |« u\ ' 

Also we use the oscillation and the total variation of the path: 

IHU,[ s ,t] = jaax \w(u) - w{v)\, (2.9) 

s<u<v<t 
N 

\\ w \\[s,t] = sup V \w(t k ) - w(t k -i)\, (2.10) 

A ti 

where A = {s = to < • • • < ijv = t} is a partition of the interval [s,t]. 

Lemma 2.3. Assume (A) and (B). Let us fix < 9 < 1. Then there exists a positive constant 
Ci such that 



M < (l+IMlg Ml£ , /2 (*-s)) e C2||w|l -^|klU )[M] . 



(2.11) 



Proof. The proof of this lemma immediately follows from the proof of Proposition 3.1 and 
Theorem 4.2 in [8]. □ 
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Note that the term ||ttf||^ r a t i g/ 2 can he replaced by a quantity defined by a modulus of 
continuity of w. We emphasize that we just need the continuity of w to estimate the bounded 
variation norm of <fi. Also we note that this estimate is not sharp in the sense that the quantity 
on the RHS does not depend on the starting point x although ||^»||[ S) t] does. If w is a continuous 
bounded variation path, we can prove the following estimate. This estimate is used to prove the 
exponential integrability of Y in the proof of Lemma 14.51 



Lemma 2.4. Assume condition (A) and the existence of the solution £ to the Skorohod problem 
for a continuous bounded variation path w. Then the total variation of the solution £ has the 
estimate: 

||£|| M <2(v / 2 + l)|k|| M (2.12) 

Proof. We write 

u(s,t) = \\w\\ M , Vo(s,t) = m)-Z(s)\, r)(s,t) = U\\ [s>t] . (2.13) 
By Lemma 2.3 (ii) in [8], 

m) - Z(s)\ 2 < \w(t) - w(s)\ 2 + - f Vo (s, u) 2 d\<P\ u + 2 /* (w(t) - w{u),d<p(u)) . (2.14) 



Noting 



ro 



— / r] (s,u) 2 d\(t>\ u < — [ r](s,u) 2 d u r](s,u)+ I rj(s,u) 7 d u u(s,u) 

< 1 Qt7(s, tf + V (s, t) 2 u(,.. I))-. k(s. I ) (2. to) 



and 



t 

2 / (w(t) -w(u),d<p(u)) 



< 2 



f 

u(u,t)d u uj(s,u) + / uj(u,t)d u r](s : u) 

ft 

(uj(s,t) — ui(s,u))d u uj(s,u) + / oj(s,t)d u r)(s,u) 



< 2 

= t) 2 + 2w(a, t)r/(s, t), (2.16) 

we obtain 

7/o (s, t) 2 < 2w(s, t) 2 + 2w(s, i)r/(s, i) + k(s, t) (2.17) 

and 

2??o(s,i) 2 < r/ (s,t) 2 + r/(s,t) 2 < w(s,t) 2 + t) + r/(s, i)) 2 + fc(a,t). (2.18) 
Therefore we have 

\/27/o(s, i) < 2w(a, i) + r/(s, t) + y/k(s,t). (2.19) 

Note that 

n 

V(s,t) = lim V]7/o(ti_i,ti), 

where A is a partition s = to < ■ ■ ■ < t n = t and |A| = snp^U — ij_i). Using the additivity, 

^(•M) = YA=iV{ti-i,ti) and w(s,i) = £)?=i u (ti-i, U) for any partition s = i < • • • < t n = t, 

we get ^/2i](s,t) < 2w(s,t) + T)(s,t) which proves the desired inequality. □ 
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Remark 2.5. Under the admissibility of the domain, Lions and Sznitman proved that ||0||r Si t] < 
HHI[«,t] which implies < 2||u;||[ S) t]. They use regularity property of the distance function 

from the boundary dD. So we may need some regularity condition on the boundary to prove 
such a stronger estimate. We note that there is a study of the regularity of the distance function, 
e.g., [7j. However, the estimate (|2.12|) is enough for our purposes. 

Let us recall the existence of strong solution and the uniqueness which is due to [HJ [5j E] . 
Let (f2, J-, P) be a complete probability space and Ft be the right-continuous filtration with the 
property that Tt contains all null set of ($7, F, P). Let B = B(t) be an J-i-Brownian motion on 
E n . Let a G C{R d -> R n <g> R d ), b G C(R d ->■ R d ) be continuous mappings. We consider an SDE 
with reflecting boundary condition on D: 

X(t)=x+ [ a(X(s))dB(s) + f b(X(s))ds + $(t), (2.20) 
Jo Jo 

where x E D. We denote this SDE by SDE(cr, b) simply. A pair of JVadapted continuous 
processes (X(t),$(t)) is called a solution to (|2.20p if the following holds. Let 

Y{t) = x+ [ a(X{s))dB(s) + I b{X(s))ds (2.21) 
Jo Jo 

Then (X(-,lj),<&(-,u)) is a solution of the Skorohod problem associated with Y(-,co) for almost 
all ijj £ fl. The following result is due to [8]. 

Theorem 2.6. Assume D satisfies condition (A) and (B) and a and b are bounded and global 
Lipschitz maps. Then there exists a unique strong solution to f|2.20|) . 

Here we note the following. This follows from the Garsia-Rodemich-Rumsey's estimate. 

Lemma 2.7. Let F = F(t,uj) be a Revalued continuous process with the property that for all 
P > 1 

E[\F(t) -F(s)\ 2p ] < C p \t- s\ p 0<s<t<T. (2.22) 

Then for all < 9 < 1 and p > 1 there exist constants C' p e which depends only on C p and 6 
such that 

^[II-^II^JO.T],^] - C 'p,6- ( 2 - 23 ) 

Since maxo<t<r \ Y(t, co)\ is exponential integrable function, using Lemma [2 .31 and Lemma f2.7l 
and Burkholder-Davis-Gundy's inequality, we immediately obtain the following estimate. 

Lemma 2.8. Assume the same assumptions as in Theorem 12.61 Let p > 1. There exists a 
positive constant C p such that 

E[\\X\\X [s , t] ]<C P \t-s\ p , (2.24) 

E[\*\\*Z t] ] < C P \t - s\ p . (2.25) 
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From now on, we always assume that a belongs to C£ and b belongs to C\. Now, we are 
going to explain our main theorem. Let N € N. Let X (t) be the solution to the reflecting 
ODE: 

X N (t) = x+ f a{X N (s))dB N (s) + I b(X N (s))ds + $ N (t), (2.26) 
Jo Jo 

where 

B N (t) = B{ttx) + ~ d) < t < t£, (2.27) 

kT 

A N B k = B(t%)-B(tti), A N = T/N, = — . (2.28) 

We already explained the existence of the strong solution to a reflecting SDE driven by a 
Brownian motion. The definition of the solution to the above equation is similar to reflecting 
SDE. However it is not trivial to see the existence of the solution. Actually, the solution to this 
equation exists uniquely by Proposition 14.11 in Section HI The following is our main theorem. In 
this paper, we do not intend to obtain the best order. The order given below is probably far 
from best. 

Theorem 2.9. Assume (A), (B) and (C). Let X be the solution to SDE(<r, b), where b = 
b + ^tr(Da)(a). Let < 9 < 1. For any p > 1, there exists a positive constant C P: x,e such that 
for all iV G N, 



E 



max \X N (t) -X(t)\ 2p 

0<t<T 



< C p>Tl eA% 6 . (2.29) 



As we noted, although this estimate may not be good, by this result and Borel-Cantelli 
lemma, we can conclude 

lim max \X l (t) - X(t) \= almost surely. (2.30) 

In order to prove this theorem, we need the Euler approximation of the solution. We explain 
the Euler approximation in the next Section. 



3 Euler approximation 

In this section, we consider the Euler approximation Xg of X. For < k < N, set = kT/N. 
Let us define X^(t) (0 < t < T) as the solution to the Skorohod problem inductively which is 
given by Xf (0) = x G R d and 

X% (t) = Xg (tti) + °{x% (tti))m) ~ S(t£U)) + K*% (*))(< " 

+ Ci<*<^- (3-1) 
In other words, X% satisfies 

X$(t) = X+ f t G(X»(7T N (s)))dB(s)+ fbiX^TT^Ws + ^it), (3.2) 

Jo Jo 
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where vrjv(t) = maxjifc | tk < t}. Define 

F#(t) = x + /" (7Ti>r( S )))dB( S ) + /" 6pf# (^(s)))^. (3.3) 

Then by the definition of the solution of the SDE, it holds that 

r(Y#)(i). (3.4) 

We prove 

Theorem 3.1. Assume (A), (B) and (C). Then for any p > 1, there exists C p > such that 



E 



max \X%{t) -X{t)\ 2p 

0<t<T 



< C p A p N . (3.5) 



To prove this theorem, we need the following lemma. 

Lemma 3.2. Assume (A) and (B). Let p > 1. There exists a positive constant C p which is 
independent of N such that 

El\\X%fV )[Sjt] ]<C p \t-s\v, (3.6) 



E 



(3.7) 



Proof. It suffices to prove (|3.7p . This follows from exponential integrability of H^b ||cc,[o,Tb 
Lemma 12.31 Lemma 12.71 and Burkholder-Davis-Gundy's inequality. □ 

Proof of Theorem \3.1[ The following proof is a modification of that of Lemma 3.1 in [5j. Note 
that we need just Lipschitz continuity of a and b and their boundedness in the proof below. It 
suffices to prove the case where p > 2. Define 

Z N (t)=X%(t)-X(t), (3.8) 

M<) = e-^ (/(Jr - ( * )+/W ' )))) , (3-9) 
k N (t) = fi N (t)\Z N (t)\ 2 . (3.10) 

Then we have 
dkM{t) 

= MJV (t)|2 (a(X# (n N (t))) - a(X(t)))dB(t)) 

+ 2 (Z N (t),b(X% (*N(t))) - b(X(t))) dt 

+ tr ((W)(Xf (7rjv(t)))) eft + tr ((W)(X(t))) di 

- tr ((V(X(i))a(X£ (**(*))))) - tr (^))V W)))) tft J 

+ 2 M7V (t) (Z*(i),d*£(t)-d$(t)) 
2fi N (t) 



1 



Z N (t)\ 2 {{(Df)(X%(t)),d$%(t)) +((Df)(X(t)),dm)} 



Z N, 



r (i)| 2 {((Af)(X£ {t))M*% (Mt)))dB(t)) +((Df)(X(t)),a(X(t))dB(t))} 
+ R N (t)dt, (3.11) 
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where 



R(t) = {(Df)(X% (t)),<x(x£ (^(t)))* - (*>#)))) (Z*(t))) eft 

4/Xjy(t) 

7 

2/ ' v( V*(*)lM((Af)(*£(^ 



+ 
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7 



(t)| 2 {trp 2 /)(^f(t))(^(Xf(7r JV (t))., C T(X#(7r i v(t))))) 



+ 



+ ti(D 2 f)(X(t))(a(X(t).,a(X(t)).))}dt 



||( J D/)(X^(i))(a(Xf(^(t)))) + ( J D/)(X( i ))(a(X(t)))f|Z-( i )|^t (3.12) 



2| ryN t + \\2. 



T 



Note that by condition (C), 

(X%(t)-X(t),d^(t)-d1>(t)) 

- ~ |Xf (t) - X(t)\ 2 { (p/)(Xf (t)), d^(t)) + ((2?/)(X(t)), (»(*))} < (3.13) 



and sup < t < T £[|Xf (t) - X% '{ir N (t))\ p ] < CA p / 2 . As for the first term on the RHS of (pTTI]) . 
using Bukholder-Davis-Gundy's inequality, we get 



E 



sup 

0<t<T 



< CE 



+ CE 



< C T E 



m (t) (Xf (t) - X(t), a(X» (Tr(t))) - a(X(t))dB(t)) 



jf \X%(t)-X(t)\Ut) 
£ |Xf (t) - X(t)| 2 |Xf (™(t)) - X# (i)| 2 di 



p/2 



k N {t) p dt 



+ C T E 



\X%{it N {t))-X%{t)\*& 



<C T [ E[k N {t) p ]dt + C T A P . 
Jo 

We can estimate the other terms similarly and we obtain 



(3.14) 



E 



sup /cjv(i) P 
0<t<T 









sup kjy(s) p 


/o 


0<s<t 



dt. 



By the Gronwall inequality, this implies the desired estimate. 

4 Proof of main theorem 

First, we prove the existence and uniqueness of the solution to reflecting ODE. 



(3.15) 
□ 



Proposition 4.1. Assume the conditions (A) and (B) hold. Let w = w(t) be a continuous 
bounded variation path on W 1 . Then there exists a unique continuous bounded variation path 
x = x{t) on R d satisfying the reflecting ODE: 

x(t) = x+ ( a(x(t))dw(t) + I b(x(t))dt + <Z>(t). (4.1) 
Jo Jo 

Proof. The following proof is a modification of the proof of Theorem 5.1 in [BJ. Note that the 
boundedness and the continuity of a and b are sufficient for the existence of the solutions. Let 
us consider the partition of [0, T] by tu = kT/N. Let x N be the Euler approximation of the 
solution, that is, let us define x N as the solution of the Skorohod problem with x N (0) = x: 

x N (t) = x N (t»_i) + *(* N (t£Li))M*) - f (tf-i)) 

+ b(x N (tt 1 ))(t - Ci) + - <^-i) <t<t». (4.2) 

Let 

y N (t)=x+ f a{x N (Tr N (s))dw(s) + f b(x N {^ N {s)))ds (4.3) 

Then {y N } is a family of uniformly bounded equicontinuous paths defined on [0, T] with values in 
W 1 . Therefore by the Arzela-Ascoli theorem, there exists a subsequence {y Nk } which converges 
in the uniform convergence topology. We denote the limit by y°°. Then by the continuity of the 
Skorohod map in Theorem E21 x Nk {= T(y Nk )), $ N k(= L(y N ") also converges to a continuous 
paths, say, <&°°, in uniform convergence topology. Clearly, the pair (x 00 ,^ 00 ) is a solution 
of a Skorohod problem associated with y°°. Taking the limit N/. — > oo in (|4.3p . we have 



y°°(t) = x + [ a(x°°(s))dw(s)+ [ b(x°°(s)))ds. (4.4) 
Jo Jo 

This shows that $°°) is a solution of the reflecting ODE. We can check the uniqueness in a 
similar manner to Theorem 5.1 in [8]. Note that the boundedness of a and b and their Lipschitz 
continuity are sufficient for the proof. □ 

Remark 4.2. We may prove the existence of the solution of reflecting ODE when the driving 
path is just p- variation path, where 1 < p < 2 using Davie's argument Q]. We will study this 
problem hopefully together with more general rough differential equation which corresponding 
to p > 2 in future's paper. 

By the proposition above, there exists a unique solution to (|2.26p . From now on, for sim- 
plicity, we may denote Aj^B^, An, Vu by ABf., A, if.. By the definition, it holds that 



X N (t) = X N (t k _ 1 ) + f a (X N {s))^ds+ f b(X N (s))ds (4.5) 

+ q>*(t) - <f> N (t k ^) 4-i < t < t k . (4.6) 
Clearly, X N (tk-i) is Tt k _^ -measurable. Let 

Y N (t) = x + f a(X N (s))dB N (s)+ f b(X N (s))ds. (4.7) 
io io 

Then X N = T{Y N ) and <S> N = L(Y N ). 
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Lemma 4.3. Assume (A) and (B). Fix JVgN. Let tk—i < s < t < t k . The constant C below 

is independent of t, s, k, N. 

(1) The following relations hold. 



and 



Y N (t)-Y N (t k _ 1 )= f a(X N (s))^ds + [ 



a(X N (s))^ds = a(X N (t k ^))^(t - t k -i) 



b(X N (s))ds (4.8) 

fc-i 



+ I I / l DaHS Ur)»K^ir) I 



f (f (Da)(X N (r))d<S> N (r) 



A 

= /* (t) +/*(*)+ If (t) + if (t). (4.9) 
Also set If (i) = / t *_ i 6(X JV (s))ds. Then 

|jf(t)| < C\AB k \^ {t ~ ^- l)2 , (4.10) 

\I^t)\<C\AB k \ {t -^~ l)2 , (4.11) 
i2 — *fc— i \ 2 . (* — ife-i) 2 x 



(2) We have 



|l£(t)| < C \ ( 1 + v ^ |AB»| , (4.12) 



|I 4 fc WI <C(t-t fc _i). (4.13) 



\Y N (t) -Y N (s)\<c(\AB k \ t -^ + t-s^ (4.14) 



^^ll^^cflA^I^-i + t-s). (4.15) 



[s,t] 

Proof. The proof of the equation (|4.8p and (|4.9p is a simple calculation. The estimate in (|4.14j) 
follows from f|4. 8|) . Hence the estimate (|4.15j) follows from this estimate and Lemma [2.4i By the 
boundedness of a, Da, b, we get KWh . KTTh . (I4TT3]) . Using (gHHJ), 



l^(*)l<e7||^|| fa _ 1|f] (t - tfc - A l)|Ag * 



A 

< C (|A5 fe | 2 (^^) 2 + (t "^- l)2 |Ag fc |^ • (4-16) 
This complete the proof. □ 
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Lemma 4.4. Assume (A) and (B). Let p > 1. There exists a positive constant C p which is 
independent of N such that for all < s < t < T, 

E[\\Y N \\X [sA ]<C P \t-s\P. (4.17) 

Proof. Pick two points < s < t < T. First consider the case where there exists 1 < k < N 
such that i fc _i < s <t <t k . Then by (f4TTi|) . 

max -^Wl < C(|AS fc |^+t-a). (4.18) 

s<u<v<t A 

Hence -E[||y JV || o,[s,t]ll 2p ] < ~ s ) p - If *fc-i < s < t k < t < t k+1 for some fc, noting 

W \\ao,[s,t] - \\ Y \\oa,[s,t k ] + 11^ lloo,[t fc ,i]> (4.19) 

we can reduce this case to the first one. We consider the other cases. Let us choose 1 < I < 
m — 1 < N such that ti_i < s < t[ < t m -i < t < t m . Then 

Y N (t)-Y N {s) 

4 r m-1 ^| 

= E 4(*o - 4oo + E - ^(^-i)) + cw - c(wi) (4.20) 

n=0 I fc=«+l J 

i 



E(^(*) - j^w). 



n=0 



Note that { J^(t) | < i < T} are continuous processes and it suffices to estimate J^[|^ r s t j]. 
First let us consider the term Jq . Let M N (t) be a continuous J^-martingale such that 

M iV (t)= [* o-(X N (ir N (s)))dB(s). (4.21) 

JO 

Then is the piecewise linear approximation of M N at the times {t k } k=1 . Therefore, 

|| Jo^lkM < , max \M N (t k )-M N {t k ,)\ 

l — l<k,k'<m 

<2 max \M N (t k ) - M N (ti)\ 

l—l<k<m 

<2 max \M N (r) - M N (ti)\. (4.22) 

i;-i<»*<tm 

Using Doob's inequality, we get 

£[H^l£ M ] < C p (t m - U.^ < $>C p (t - sf. (4.23) 
Next we consider the term . By the estimate (|4.12p . we have 

ll^llocM < C ^E \^ B k\ 2 + A • |A£ fc |^ . (4.24) 
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Note that 

{AB k }f =l = VX{Z k }Z =l inlaw, (4.25) 
where {6c}fcL; are i.i.d. random variables whose common distribution is N(0, 1). Therefore 



E 



N\\2p 

oo, [s,t] 



< C p A 2p E 



< C p A 2p E 



' m \ 2 P 
\k=l / 



+ C p A 4p (m -l + l) 



+ C p A 4p (m -l + l) 



2p 



< C p A 2p + Cp A4p (m -l + l) 



Kk=l 



Cn 



m 



Z + 1 \ 2p 



T < CJt 



\2P 



N J 

We can estimate other terms in a similar way and we complete the proof. 



(4.26) 



□ 



Lemma 4.5. Assume (A) and (B). Let p > 1. There exists a positive number C p which is 
independent of iV such that for all < s < t < T, 



mx N C, M ]<Cp\t-'F, 

E[\<S> N \\ 2p t] ]<Cp\t-s\ p . 



(4.27) 
(4.28) 



Proof. It suffices to prove (|4.28p . By checking the exponential integrability of \\Y ||oo,[o,t]j we 
can prove this by using the fact & N = L(Y N ), Lemma 12.31 Lemma 14.41 and Lemma 12.71 Note 
that 



E 



N 



exp 



\k=l 



9 T\ ~ nN / 2 

1-^-) -f" r " as.V-x. 



(4.29) 



This and the estimates in Lemma 14.31 (1) and ()4.22p implies the exponential integrability of 



\Y 



loo,[0,T]- 



□ 



The following is a key lemma for the proof of LP convergence of Wong-Zakai approximation. 

Lemma 4.6. Assume (A), (B) and (C). Let be the Euler approximation to SDE(cr, b), 
where b = b + |tr (Da)(a). Then for any < 9 < 1, there exists a positive constant Cq such 
that for all N, 



6/2 



sup E[\X-(t-)-X-(t^)\ 2 \<C e .A N 

0<k<N 



(4.30) 



In the proof of this lemma, the integrals which contains J^-semimartingales and non-adapted 
bounded variation processes, e.g. Wong-Zakai approximation Xjy(t) appear. Hence we need the 
following definition of the integrals. 
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Lemma 4.7. Let X(t),Y(t) be J-t-continuous semimartingales and A(t) be bounded variation 
continuous process. Suppose that sup 0<t<T {\X (t)\ + |Y(i)|} + |-*4( - )l[o,T] £ L p for all p. Define 

/ X(s)A( S )dY(s) = lim ^X(tf_ 1 )A(^_ 1 )(F(^)-y(tf_ 1 )), (4.31) 



fc=i 



(XA,Y) t = lim £ - (X^O^f-l))) (^) " n<£L0), 

N— too *■ — ' 

k=l 

(4.32) 

where = tk/N. These converge in probability and it holds that 

/ X(s)A(s)dY(s) = [ A(s)dZ(s) = A(t)Z(t) - [ Z{s)dA{s), (4.33) 
Jo Jo Jo 

(XA,Y) t = [ A(s)d(X,Y) s , (4.34) 
Jo 

where Z(s) = /„* X(s)dY(s) is usual Ito integral and the RHS of ()4.33p is Riemann-Stieltjes 
integral. 

Let us consider a set of stochastic processes § which consists of a finite sum of product process 
Y(t)A(t), where Y{t) is a Tf continuous semimartingale and A(t) is a continuous bounded vari- 
ation process which is not necessarily JVadapted with the property sup <£<T 1^(^)1 + H-^H [o,T] £ 
C] p >iL p . Then this class is stable under the stochastic integral in the sense of lemma above. In 
the calculation below, we use the integrals of stochastic processes in this sense. Moreover the 
following chain rule holds. 

Lemma 4.8. Let Y, Z G §. Then 

Y{t)Z{t) = r(0)Z(0) + f Y(s)dZ(s) + I Z{s)dY(s) + (Y, Z) t , (4.35) 

Jo Jo 

where (Y, Z)t is defined similarly to Lemma 14.71 

The above two lemmas are proved by a standard argument and we omit the proof. The 
proof of Lemma 14.61 is long. In that proof, we use the following two lemmas several times. 
First lemma is an estimate of the expectation of stochastic integrals above. To this end, we 
introduce a family of iterated integrals. Let S be a set of stochastic processes which consists of 
the processes g(Y(t)) where g is a C 1 function with values in R with bounded derivative and 

Y = X N ,X%,B,B N ,$ N (t),$% (t). (4.36) 

We define a set Si of two parameter processes / = f(s,t) (0 < s < t < T) inductively. Let 
So = {!}• The set Si (i > 1) consists of finite sums of 

rt 



llfk(s,t), [ g(s,u)df (u), (4.37) 
where f k £ S ik J2i=i i k = i h > 1 and f £S,g £ 



k=l 
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Lemma 4.9. Let t^_ l < s < t < . For any / € Si (i G N), there exists C p > which is 
independent of N, k such that 

|| max f(u,v)\\ LP <C(p)(t-sr /2 . (4.38) 

s<M<?;<i 

Proof. In this proof, we say that / G Uj>o5i is adapted when the following holds. The definition 
is given inductively by 

(i) 1 € 5o is adapted, 

(ii) Let / be finite linear sums of processes in (|4.37|) . Then / is adapted if all (1 < A; < j) 
and g are adapted and f = g{Y(t)), where Y = X$,B,$% and g is a C 1 function with 
bounded derivative. 

It is easy to check that the set Si is equal to the set of finite sums of two parameter processes 

(jl [ u)dA k {u)^ ■ h(s, t). (4.39) 

Here A k is a bounded variation process in 5 and E Si k . When p = 0, we set this term as 
1. Also h £ Sj is adapted and satisfies Ylk=i^ k + 1) + J = ^- This easily follows from the 
induction on i. Using this result, Lemma 13.21 and Lemma 14.51 we can complete the proof of the 
desired result by an induction oni. □ 

Proof of Lemma 14.61 We write 

Z N (t) = X%(t)-X N (t), (4.40) 

PN (i) = e-^ X ^ t)+ ^ XNm \ (4.41) 
m N (t) = PN (t)\Z N (t)\ 2 . (4.42) 
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Let tk—i <t<t k . By Lemma 14.81 

dm at (t) 

AB h 



PN (t)l2 \Z N (t),a(X§(t k ^))dB(t) - a (X N (t))^dt 

+ 2 (z N {t)MX%{t k ^)) - b{X N (t))) dt + tr [(Va)(Xf (t k ^))] dtj 
+ 2 PN (t) {Z N (t),d^(t)-d1> N (t)) 

- ^ \Z N (t)\ 2 {((Df)(X% (t)),d*%(t)) + ((Df)(X N (t)U* N (t))} 
2/ ' V(/) \Z N (t)\ 2 {((Df)(X%(t)),a(X»(t k ^))dB(t)) 
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+ ((D/)(^(t)) l ^(t))^)} 
^1 {{Df){X% (t)), (tjfc-i))*) <r(*f (t fc _i))e<) tft 

2pjv(t) 



7 
7 



Z w (t)| 2 ^(D/)^^))^^^-!))) + ((£>/)(X^ (t)), 6^ (*))) eft 
Z"(t)| 2 tr(£> 2 /)(*£ (*)) (t*-i)-,a(Xf (t fc _!))-] (it 



7 2 

where {e^} is a c.o.n.s of M n . After integrating both sides from t k -\ to t k , we see that the sum of 
the integral of the second term and the third term on the RHS is non-positive by the condition 
(C), Therefore 

6 

m N (t k ) < m N (tk-i) + ( 4 - 44 ) 

k=l 

where 

h = jf* PN {t)U {z N {t),a{X^{t k ^))dB{t) - a (X N (t))^dt 

+ 2 (z N (t),b(X% (t fc _!)) - b(X N (t))) dt + tr ((Va)(Xf (t fc _ a ))) dtj (4.45) 

Yi( D f)( X E (t)),*(XE (ife-i)K) (^(i),^ (t fc _i))e<)di (4.46) 



t h -i 7 

'' 2 

tk-i 7 



J: 



3 = - -mN(t){((Df)(X%(t)),a(X%(t k ^))dB(t)) 



+ [(Df)(X N (t)),a(X N (t))^dt)\ ( 1. IT) 



0} 
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h = ~ f k -m^(t)(((Z)/)(Xf(t)),6(Xf(t fc _ 1 )))dt+(p/)(X^(t)),6(X^(t)))dt) (4.48) 



/*" ^^trp 2 /)(^f (*)) [*(*f (tfc-i)-.a(Xf (i fe -l))-] dt (4.49) 



1 6 



^^KZ?/)^^))^^-!)))! 3 *. (4.50) 

Let a/j = E[m]y(tk)]- We prove that there exists a positive constant C and < 9 < 1 which is 
independent of N and a non-negative sequence such that 

a fe < (l + a fe _i + 6 fc 1 < fc < iV (4.51) 

N 6/2 



k=l 

Then we get 



«fc < ^1 + ^ + (l + jpj h-i + h 

fc / T x 0/2 

<e c 26i<C - (4.53) 
i=l ^ ' 

which is the desired estimate. As for //% (A; = 4, 5, 6), by Lemma 14.91 

\E[I k ]\ < ^ flfc _ 1 + A 3 / 2 . (4.54) 

So our task is to estimate Ii, 12,1s- To estimate Ji, applying Lemma 14.81 to \Z N (t)\ 2 , we have 

h = Ji + J 2 + J 3 + J 4 , (4.55) 

where 

Ji = 2 jT* ^(t fc -i)| (z N (t),v(XE (tk-i))dB(t) - a(X N (t k ^))^dt^ 

+ itr(W)(Xf(t fc _i))cft| (4.56) 

(4.57) 
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J 3 = 2^ ( PN (t) - pjv(* fc _i)) (V(t),cT(X# (t fe _ 1 ))d J B(t)- ( r(X Ar (t jfc _ 1 ))^^ , (4.58) 
J 4 = / * (p JV (t)-/9iv(^-i))tr(W)(Xf(t fe _ 1 ))^. (4.59) 



We estimate each terms J«. First we estimate J\. Let 



Jx = 2 jT* p N {t k ^)^Z N {t) - Z N (t k _ l ),a{X^{t k _ l ))dB{t) - a (X N (t k _ 1 ))^d?j 

+ ^tr(W)(Xf(t fe _ 1 ))^|. (4.60) 
Then i?[Ji — Ji] = 0. So it suffices to estimate the expectation of J\. We rewrite 

4 

Ji = ^Ji,k, (4-61) 



fc=i 

where 



rtk ( / Ag 

Ji,i = 2 jf ^(tfc-i)j (t*_i))(5(t) - B(t t -i)) - ^(t*.!))-^^ " **-i)> 

AD v 

a(X%(t k ^))dB(t) - a(X N (t k ^))^dt) 

+ hr( t aa)(X%(t k _ 1 ))dt\. (4.62) 



Ji, 2 = -2 



*fc- 



P7v(tfc-i)|(^ (<x(X%)) - a (^(tfc-x))) ^ds, 



(t fc _!))dB(t) - a (X N (t k _ 1 ))^dt) | (4.63) 
Ji,3 = 2f tk p N (t k -i)\ (b(X%(t k _ 1 )))(t-t k - 1 )- f b(X N (s))ds, 



<r(X#(t* ; ;.;u//,'.;/i rr.;.V v .;/, ^ (LGl) 



J M = 2 J*" ((<&£ (t) - ^(t fc -i)) - - ^(tfc-l)), 



(* fc _i))dB(i) - a{X N {t k ^))-^dt 
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>} 



(4.65) 



By a simple calculation, 

E[J 1A ] = E [||a(Xf (**_!)) - a(X N (t k ^))\\ 2 H _ s ] (t k - t k -i) < Ca k A. (4.66) 
By Lemma |4.9| we have 

E[\Jx,2\] < CA 3/2 . (4.67) 
It is easy to see that E[\ J13I] < CA 3 / 2 . By integrating by parts 

\E[J 1A ]\ < CE [{\\$%\\ [th _ 1>tk] + [IBIU,^^] =: c k . (4.68) 

We have 



N 



k=l 

We estimate J2 

"ft 



J> < Cll llM + ll^ll[o,T]) ll^llmaxHSIU,^^!!^ < CA e ' 2 . (4.69) 



Jt k -i v ' 

+ 2 f p N (t)(z N (t),b(X N (t k - 1 ))-b(X N (t k - 1 )) 

Jtk-1 v 

-(a(X N (t))-a(X N (t k ^))) ^)dt 

= J 2 ,i + J2,2- (4.70) 

By rewriting = Z N (t k -i) + - Z N (t k ^) and using Lemma 14.91 and the Schwarz 

inequality, we get 

E[\J 2 ,i\] <C(a k A + A^ 2 ). (4.71) 

We consider J 2 ,2- Let 

J 2 , 2 = 2 / * PN (t k „ 1 )(z N (t k „i)MX N (t k -i))-b(X N (t k „ 1 )) 

- (a(X N (t)) - a(X N (t k ^))) ^)dt. (4.72) 

By Lemma |4.9| we have 

\E[J 2 ,2-J2,2]\<CA 3 / 2 . (4.73) 

Noting 

a(X N (t))-a(X N (t k ^)) 

(Da)(X N (s)) (a(X N (s))^J ds + £ (Da)(X N (s))(b(X N (s)))ds 

+ f {Da)(X N \s)){d$ N \s)) (4.74) 
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and for any £, 



E 

we obtain 
Since 



{t-t k ^){Dam{a^)^j[^\dt 



-tr(Da)(0(a(0) = (b-b)(0, (4.75) 



\E[J 2;2 ]\ < CE [||^ JV ||[t fc _ 1 ,t fe ]||-B||oo,[t fc _ 1 ,t fe ]] +CA 3 / 2 . 



(4.76) 



N 



< CA 9 ' 2 



k=i 



||^|| [0)T] max|| J B|| OO)te _ ljtfe] 



(4.77) 



we obtain the desired estimate for J 2>2 . By Lemma 14.91 we have £"[1 J4I] < CA 3 / 2 . We estimate 
J3 together with I 2 and a term which is defined below. We estimate I3. We rewrite 



^3 — -^3,1 + h,2 + ^3,3 + ^3,4 + -^3,5 + -^3,6 + h,7, 



(4.78) 



where 



4 f tk 



7 Jt k _, 



PN(t) 



£ (z N (s) - Z N (t k ^),a(X^ (t fc _i))dB( 8 ) - a(X N (s))^ds^J 



AB h 



x { ((Df)(Xg{t))MXS(th-i))dB(t)) + ( (Df)(X N (t)),a(X N (t))^dt) }, (4.79) 



h,2 = — / PN(t) / ' - A 
7 Jt k _ 1 Jt k -i 



Z N (t k ^), a(X% (t fc _i))dB(a) - a(X%))=pd S 



A5 fc 
A 



x I ((Df)(X^(t)), a(X^ (tk-i))dB(t)) + ^(Df)(X N (t)),a(X N (t))^^-dt^j 1, (4.80) 



j 3> 3=-- r Piv(t) r (^(^^(^(t^oj-K^w))^ 

x j((Z>/)(X# m*(X§ (**_i))dB(t)) + ^{Df)(X N (i)),a (X N (t) ) ) (4.81) 



I 3) 4 = r PJVW / (^(S) " Z N (t k -!),d^(8) - d$ N ( S )) ds 

7 Jt k _\ Jt k -i 



x I ((Df)(X^ (t)),a(X^ (t k _i))dB(t)) + ^ (Df)(X(t)), a(X (t))—^-dt 



AB k 



(4.82) 
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7 



(i)),<r(X£ (t*_i))dB(t)) + [, (4.83) 



/ 3 , 6 = --|^(t fc -i)| 2 T PN(t)\((Df)(X^(t)),a(X^(t k ^))dB(t)) 

7 



(4.84) 



4,7 = -- -t fc _i)tr((W)(X# (t fc _0)) {((^(Xf^))^^!^!))^^)) 

7 -/tfc-1 



AS, 



+ ((JV)(X"(t)),a(jr"(t))-^dt 



As for ^3,1, ^3,3, -^3,4, ^3,7, by Lemma I4~9l it is easy to see 

\E[I 3>k }\ <CA 3 / 2 . 



(4.85) 



(4.86) 



Also similarly, i^/3^1] < Aajt-i- Using Lemma 14.91 we have there exists non-negative random 
variable I' 3 5 such that E[I 3 5 ] < CA 3 / 2 and 



|/ 3l5 | <C\Z N {t k ^)\G k + l' 



(4.87) 



where 



Gfc = max 

\ tk_i<t<t k 



{(Df)(X%(s)),a(X%(t k ^))dB(s)) 



+ 1151 



°°,[tfc-i,tfc] 



We obtain 
Since 



£[|/ 3 , 5 |] < C(a k ^A + £ [G fc ]) + CA 3 / 2 . 



(4.88) 



(4.89) 



fc=i 



$£||[0,T] + II<^II[0,T]) 

(( J D/)(X#( S )), f r(X#(t fc _ 1 ))di?( S )) 



x max max 

k \tk-i<t<t k 



tk-l 



+ 1151 



oo>[tfe-l,*fc] 



< CA / 2 , 



(4.90) 
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we obtain the desired estimate for I3 5. Finally we estimate I2 + J3 + ^3,2- First we rewrite J3. 
Note that 



p N (t) - p N (t k -i) = --PN(tk-i)\ ((Df)(X% (t h -i)),*(X% (t fc _i))(B(t) 



+ ( (Df)(X N (tk-i)),cr(X N (tk-i))^^{t - tk-x) 



- -Piv(ife-i) I ((£>/) (Xf &_!)), 6(Xf (t fc _i))(t - t fc _i) 
+ ((£>/)(**(**_!)), b(X N (t k -i))(t - i fe -i)) I 

" ~ P iv(i fe -l){((IV)(*# (t*_i)),<l>£ (t) - *£(**-!)) 

+ [(Df)(X N (tk-x)), ® N (t) - ^(tfc-i))! +p(t k - 1 ,t). (4.91) 

Here p £ 1S2. Hence we can neglect the term p to estimate J3 by Lemma 14.91 Also we can 
estimate the terms containing <frg,& N in a similar way to Ji,4, ^2,21 -^3,5- We can estimate the 
term containing 6, 6 by Lemma 14.91 Consequently, we can replace the term J3 by J3: 



4 



k 



■h - ~ jf pAr(i*-i)j (tfe-i)),<r(Xf (t k -x))(B{t) -B(t*_i)) 



+ ( ( J D/)(X JV (t fc _ 1 )), ( 7(X JV (i fe _ 1 ))^(t-t fe _ 1 ) 

x ( ^(tfc.O.a^f (t fc _!))dS(t) -cT(X N (t k -x))^dt ) . (4.92) 

Also, similarly, we can replace 13^ by ^3^: 
4 

^3,2 = / PN(tk-x) 

x (V^-x), a(x£ (t fc _ x ))(B(t) - fl(* fc -i)) - a (X N (t k ^))^(t - t fc _x) 
x J (0D/)(Xf (tfc.!)), a(Xf (t k _x))dB(t)) + ((Df)(X N (t k _ 1 )),a(X N (t k ^ 1 ))^dt 



(4.93) 
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By a simple calculation, we have 



E 



J 3 I -?Vi 
OA 

— PN (t k -i) (( D f)( X E (tfc-i)), (7(Xf (t fc _i))ei) (^(tfc-i), a(X JV (t fe _ 1 ))e i ) 



7 



2A 



— PN{tk-i) ^2 {(D f)(X N (t k _i)) , a(X N (t k _i))ei) {Z N (t k ^), a(Xf (t fc _i))ei) 

OA 

+ —PN(tk-i) ^2 {{D f)(X N (t k -i)) , a(X N (t k -i))ei) (Z N (t k ^), a(X N (t k ^)) ei ) , (4.94) 



^3,2 | ^i fc _i 

OA 

^ — Piv(t fc -i) (^(t fc -i),«r(Xf (tfc.xJJei) [(Df)(X N (t k -i)), a(X N (tk-i))ei) 

OA 

+ — ^(t fc -i) (^(tfc-i), a^fa-i))*) ((Df)(X% (t k ^)),a(X% (t fc _ 1 ))e l ) 
OA 

+ — PN (t k -i) (Z N (t k ^),a(X N (t k ^)) ei ) ((Df)(X N (t k ^)),a(X N (t k ^)) ei ) . (4.95) 



Hence 



^3 + 4,2 

Ofc_iO(A) 



4A „ 
+ — £ 

7 



Consequently, 



p N (t k - 1 )^2(Z N (t k . 1 ),a(^ N (tk-i))e i ) {(Df)(X N (t k ^)),a(X N (t k ^)) ei ) 



E 



h + h + h:. 



< Ca fe _iA + A 3 / 2 . 



This completes the proof. 



(4.96) 

(4.97) 
□ 



Remark 4.10. Note that some parts in the above estimate for a k = E[m n (t k )] is crude. In the 
case where dD = 0, that is, D = R d , the local time term <1>^ = & N vanish. In this case, the 
estimate 



a k < afc_iA + CA 2 



(4.98) 



might be true. The bad term A e / 2 essentially comes from the estimates on local time terms if 
dD + 0. 

Lemma 4.11. Assume the same assumptions in Lemma 14.61 and consider the same SDE. Let 
< 8 < 1. Then there exists a positive constant C p ,T0 such that 



E 



max \X N (t)-XS(t)\ 2p 

0<t<T ' W h Wl 



< r a 9 / 6 



(4.99) 
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Proof. Let Nq be a natural number and choose a sufficiently large natural number N. Pick 
partition points {sfc}^ x C {t^}^ =1 such that 



l*fc — 



< 



N' 



(4.100) 



Let <t<t%°. Then 



\X N (t) - X% (t)\ < \X N (t) - X N (t^)\ + \X N (t»° ) - X N (s k )\ + \X N (s k ) - X" (s k ) 



and 



+ \X% (s k ) - X% (C)| + \X% (C) " (*)l 



(,N 



(4.101) 



max -X£(t)| < 2 



0<t<T 



max 

0<s<t<T,|t-s|<T/iVo 



i w (t)-r( s 



Therefore 
E 



max iX^m - A^(i)| 2p 

0<i<T 



< 



+ 2 max \XS (t) - XU is) 

0<s<t<T,\t-s\<T/N o 
N 

+ J2\X N (s k )-X%(s k )\. 

k=l 



(4.102) 



N 's k W 



+ N 2 n 



/ rp \ 0/2 



(4.103) 



Here we use the uniform moment estimate for A^ , A^ and Lemma 14.61 and Lemma 12.71 Hence 
setting iVo as the integer part of A rl / 6p , we obtain the desired estimate. □ 

Proof of main theorem. The proof follows from Theorem 13.11 and Lemma 14. Ill □ 
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